Epigenetic Reprogramming of Autophagy Drives Mutant IDH1 Glioma Progression and Response to Radiation.

Mutant isocitrate dehydrogenase 1 (mIDH1; IDH1 R132H ) exhibits a gain of function mutation enabling 2-hydroxyglutarate (2HG) production. 2HG inhibits DNA and histone demethylases, inducing epigenetic reprogramming and corresponding changes to the transcriptome. We previously demonstrated 2HG-mediated epigenetic reprogramming enhances DNA-damage response and confers radioresistance in mIDH1 gliomas harboring p53 and ATRX loss of function mutations. In this study, RNA-seq and ChIP-seq data revealed human and mouse mIDH1 glioma neurospheres have downregulated gene ontologies related to mitochondrial metabolism and upregulated autophagy. Further analysis revealed that the decreased mitochondrial metabolism was paralleled by a decrease in glycolysis, rendering autophagy as a source of energy in mIDH1 glioma cells. Analysis of autophagy pathways showed that mIDH1 glioma cells exhibited increased expression of pULK1-S555 and enhanced LC3 I/II conversion, indicating augmented autophagy activity. This dependence is reflected by increased sensitivity of mIDH1 glioma cells to autophagy inhibition. Blocking autophagy selectively impairs the growth of cultured mIDH1 glioma cells but not wild-type IDH1 (wtIDH1) glioma cells. Targeting autophagy by systemic administration of synthetic protein nanoparticles packaged with siRNA targeting Atg7 (SPNP-siRNA-Atg7) sensitized mIDH1 glioma cells to radiation-induced cell death, resulting in tumor regression, long-term survival, and immunological memory, when used in combination with IR. Our results indicate autophagy as a critical pathway for survival and maintenance of mIDH1 glioma cells, a strategy that has significant potential for future clinical translation. One Sentence Summary: The inhibition of autophagy sensitizes mIDH1 glioma cells to radiation, thus creating a promising therapeutic strategy for mIDH1 glioma patients. Graphical abstract: Our genetically engineered mIDH1 mouse glioma model harbors IDH1 R132H in the context of ATRX and TP53 knockdown. The production of 2-HG elicited an epigenetic reprogramming associated with a disruption in mitochondrial activity and an enhancement of autophagy in mIDH1 glioma cells. Autophagy is a mechanism involved in cell homeostasis related with cell survival under energetic stress and DNA damage protection. Autophagy has been associated with radio resistance. The inhibition of autophagy thus radio sensitizes mIDH1 glioma cells and enhances survival of mIDH1 glioma-bearing mice, representing a novel therapeutic target for this glioma subtype with potential applicability in combined clinical strategies.

Gene Target Prediction of Environmental Chemicals Using Coupled Matrix-Matrix Completion.

Human exposure to toxic chemicals presents a huge health burden. Key to understanding chemical toxicity is knowledge of the molecular target(s) of the chemicals. Because a comprehensive safety assessment for all chemicals is infeasible due to limited resources, a robust computational method for discovering targets of environmental exposures is a promising direction for public health research. In this study, we implemented a novel matrix completion algorithm named coupled matrix-matrix completion (CMMC) for predicting direct and indirect exposome-target interactions, which exploits the vast amount of accumulated data regarding chemical exposures and their molecular targets. Our approach achieved an AUC of 0.89 on a benchmark data set generated using data from the Comparative Toxicogenomics Database. Our case studies with bisphenol A and its analogues, PFAS, dioxins, PCBs, and VOCs show that CMMC can be used to accurately predict molecular targets of novel chemicals without any prior bioactivity knowledge. Our results demonstrate the feasibility and promise of computationally predicting environmental chemical-target interactions to efficiently prioritize chemicals in hazard identification and risk assessment.

Translational toxicoepigenetic Meta-Analyses identify homologous gene DNA methylation reprogramming following developmental phthalate and lead exposure in mouse and human offspring.

Although toxicology uses animal models to represent real-world human health scenarios, a critical translational gap between laboratory-based studies and epidemiology remains. In this study, we aimed to understand the toxicoepigenetic effects on DNA methylation after developmental exposure to two common toxicants, the phthalate di(2-ethylhexyl) phthalate (DEHP) and the metal lead (Pb), using a translational paradigm that selected candidate genes from a mouse study and assessed them in four human birth cohorts. Data from mouse offspring developmentally exposed to DEHP, Pb, or control were used to identify genes with sex-specific sites with differential DNA methylation at postnatal day 21. Associations of human infant DNA methylation in homologous mouse genes with prenatal DEHP or Pb were examined with a meta-analysis. Differential methylation was observed on 6 cytosines (adjusted-p < 0.05) and 90 regions (adjusted-p < 0.001). This translational approach offers a unique method that can detect conserved epigenetic differences that are developmentally susceptible to environmental toxicants.

piOxi database: a web resource of germline and somatic tissue piRNAs identified by chemical oxidation.

PIWI-interacting RNAs (piRNAs) are a class of small non-coding RNAs that are highly expressed and extensively studied from the germline. piRNAs associate with PIWI proteins to maintain DNA methylation for transposon silencing and transcriptional gene regulation for genomic stability. Mature germline piRNAs have distinct characteristics including a 24- to 32-nucleotide length and a 2'-O-methylation signature at the 3' end. Although recent studies have identified piRNAs in somatic tissues, they remain poorly characterized. For example, we recently demonstrated notable expression of piRNA in the murine soma, and while overall expression was lower than that of the germline, unique characteristics suggested tissue-specific functions of this class. While currently available databases commonly use length and association with PIWI proteins to identify piRNA, few have included a chemical oxidation method that detects piRNA based on its 3' modification. This method leads to reproducible and rigorous data processing when coupled with next-generation sequencing and bioinformatics analysis. Here, we introduce piOxi DB, a user-friendly web resource that provides a comprehensive analysis of piRNA, generated exclusively through sodium periodate treatment of small RNA. The current version of piOxi DB includes 435 749 germline and 9828 somatic piRNA sequences robustly identified from M. musculus, M. fascicularis and H. sapiens. The database provides species- and tissue-specific data that are further analyzed according to chromosome location and correspondence to gene and repetitive elements. piOxi DB is an informative tool to assist broad research applications in the fields of RNA biology, cancer biology, environmental toxicology and beyond. Database URL:  https://pioxidb.dcmb.med.umich.edu/.

Salivary microbiome changes distinguish response to chemoradiotherapy in patients with oral cancer.

BACKGROUND: Oral squamous cell carcinoma (SCC) is associated with oral microbial dysbiosis. In this unique study, we compared pre- to post-treatment salivary microbiome in patients with SCC by 16S rRNA gene sequencing and examined how microbiome changes correlated with the expression of an anti-microbial protein. RESULTS: Treatment of SCC was associated with a reduction in overall bacterial richness and diversity. There were significant changes in the microbial community structure, including a decrease in the abundance of Porphyromonaceae and Prevotellaceae and an increase in Lactobacillaceae. There were also significant changes in the microbial community structure before and after treatment with chemoradiotherapy, but not with surgery alone. In patients treated with chemoradiotherapy alone, several bacterial populations were differentially abundant between responders and non-responders before and after therapy. Microbiome changes were associated with a change in the expression of DMBT1, an anti-microbial protein in human saliva. Additionally, we found that salivary DMBT1, which increases after treatment, could serve as a post-treatment salivary biomarker that links to microbial changes. Specifically, post-treatment increases in human salivary DMBT1 correlated with increased abundance of Gemella spp., Pasteurellaceae spp., Lactobacillus spp., and Oribacterium spp. This is the first longitudinal study to investigate treatment-associated changes (chemoradiotherapy and surgery) in the oral microbiome in patients with SCC along with changes in expression of an anti-microbial protein in saliva. CONCLUSIONS: The composition of the oral microbiota may predict treatment responses; salivary DMBT1 may have a role in modulating the oral microbiome in patients with SCC. After completion of treatment, 6 months after diagnosis, patients had a less diverse and less rich oral microbiome. Leptotrichia was a highly prevalent bacteria genus associated with disease. Expression of DMBT1 was higher after treatment and associated with microbiome changes, the most prominent genus being Gemella Video Abstract.

Effects of Developmental Lead and Phthalate Exposures on DNA Methylation in Adult Mouse Blood, Brain, and Liver Identifies Tissue- and Sex-Specific Changes with Implications for Genomic Imprinting.

Background: Maternal exposure to environmental chemicals can cause adverse health effects in offspring. Mounting evidence supports that these effects are influenced, at least in part, by epigenetic modifications. Objective: We examined tissue- and sex-specific changes in DNA methylation (DNAm) associated with human-relevant lead (Pb) and di(2-ethylhexyl) phthalate (DEHP) exposure during perinatal development in cerebral cortex, blood, and liver. Methods: Female mice were exposed to human relevant doses of either Pb (32ppm) via drinking water or DEHP (5 mg/kg-day) via chow for two weeks prior to mating through offspring weaning. Whole genome bisulfite sequencing (WGBS) was utilized to examine DNAm changes in offspring cortex, blood, and liver at 5 months of age. Metilene and methylSig were used to identify differentially methylated regions (DMRs). Annotatr and Chipenrich were used for genomic annotations and geneset enrichment tests of DMRs, respectively. Results: The cortex contained the majority of DMRs associated with Pb (69%) and DEHP (58%) exposure. The cortex also contained the greatest degree of overlap in DMR signatures between sexes (n = 17 and 14 DMRs with Pb and DEHP exposure, respectively) and exposure types (n = 79 and 47 DMRs in males and females, respectively). In all tissues, detected DMRs were preferentially found at genomic regions associated with gene expression regulation (e.g., CpG islands and shores, 5' UTRs, promoters, and exons). An analysis of GO terms associated with DMR-containing genes identified imprinted genes to be impacted by both Pb and DEHP exposure. Of these, Gnas and Grb10 contained DMRs across tissues, sexes, and exposures. DMRs were enriched in the imprinting control regions (ICRs) of Gnas and Grb10, with 15 and 17 ICR-located DMRs across cortex, blood, and liver in each gene, respectively. The ICRs were also the location of DMRs replicated across target and surrogate tissues, suggesting epigenetic changes these regions may be potentially viable biomarkers. Conclusions: We observed Pb- and DEHP-specific DNAm changes in cortex, blood, and liver, and the greatest degree of overlap in DMR signatures was seen between exposures followed by sex and tissue type. DNAm at imprinted control regions was altered by both Pb and DEHP, highlighting the susceptibility of genomic imprinting to these exposures during the perinatal window of development.

Clinical, morphologic and molecular heterogeneity of HPV-associated oropharyngeal cancer.

The incidence of human papillomavirus-positive (HPV+) oropharyngeal squamous cell carcinoma (OPSCC) is rising rapidly and has exceeded cervical cancer to become the most common HPV-induced cancer in developed countries. Since patients with HPV + OPSCC respond very favorably to standard aggressive treatment, the emphasis has changed to reducing treatment intensity. However, recent multi-center clinical trials failed to show non-inferiority of de-escalation strategies on a population basis, highlighting the need to select low-risk patients likely to respond to de-intensified treatments. In contrast, there is a substantial proportion of patients who develop recurrent disease despite aggressive therapy. This supports that HPV + OPSCC is not a homogeneous disease, but comprises distinct subtypes with clinical and biological variations. The overall goal for this review is to identify biomarkers for HPV + OPSCC that may be relevant for patient stratification for personalized treatment. We discuss HPV + OPSCC as a heterogeneous disease from multifaceted perspectives including clinical behavior, tumor morphology, and molecular phenotype. Molecular profiling from bulk tumors as well as single-cell sequencing data are discussed as potential driving factors of heterogeneity between tumor subgroups. Finally, we evaluate key challenges that may impede in-depth investigations of HPV + OPSCC heterogeneity and outline potential future directions, including a section on racial and ethnic differences.

Whole-genome CpG-resolution DNA Methylation Profiling of HNSCC Reveals Distinct Mechanisms of Carcinogenesis for Fine-scale HPV+ Cancer Subtypes.

DNA methylation is a vital early step in carcinogenesis. Most findings of aberrant DNA methylation in head and neck squamous cell carcinomas (HNSCC) are array based with limited coverage and resolution, and mainly explored by human papillomavirus (HPV) status, ignoring the high heterogeneity of this disease. In this study, we performed whole-genome bisulfite sequencing on a well-studied HNSCC cohort (n = 36) and investigated the methylation changes between fine-scaled HNSCC subtypes in relation to genomic instability, repetitive elements, gene expression, and key carcinogenic pathways. The previously observed hypermethylation phenotype in HPV-positive (HPV+) tumors compared with HPV-negative tumors was robustly present in the immune-strong (IMU) HPV+ subtype but absent in the highly keratinized (KRT) HPV+ subtype. Methylation levels of IMU tumors were significantly higher in repetitive elements, and methylation showed a significant correlation with genomic stability, consistent with the IMU subtype having more genomic stability and better prognosis. Expression quantitative trait methylation (cis-eQTM) analysis revealed extensive functionally-relevant differences, and differential methylation pathway analysis recapitulated gene expression pathway differences between subtypes. Consistent with their characteristics, KRT and HPV-negative tumors had high regulatory potential for multiple regulators of keratinocyte differentiation, which positively correlated with an expression-based keratinization score. Together, our findings revealed distinct mechanisms of carcinogenesis between subtypes in HPV+ HNSCC and uncovered previously ignored epigenomic differences and clinical implications, illustrating the importance of fine-scale subtype analysis in cancer. Significance: This study revealed that the previously observed hypermethylation of HPV(+) HNSCC is due solely to the IMU subtype, illustrating the importance of fine-scale subtype analysis in such a heterogeneous disease. Particularly, IMU has significantly higher methylation of transposable elements, which can be tested as a prognosis biomarker in future translational studies.

H3.3-G34R Mutation-Mediated Epigenetic Reprogramming Leads to Enhanced Efficacy of Immune Stimulatory Gene Therapy in Pediatric High-Grade Gliomas.

Pediatric high-grade gliomas (pHGGs) are diffuse and highly aggressive CNS tumors which remain incurable, with a 5-year overall survival of less than 20%. Within glioma, mutations in the genes encoding the histones H3.1 and H3.3 have been discovered to be age-restricted and specific of pHGGs. This work focuses on the study of pHGGs harboring the H3.3-G34R mutation. H3.3-G34R tumors represent the 9-15% of pHGGs, are restricted to the cerebral hemispheres, and are found predominantly in the adolescent population (median 15.0 years). We have utilized a genetically engineered immunocompetent mouse model for this subtype of pHGG generated via the Sleeping Beauty-transposon system. The analysis of H3.3-G34R genetically engineered brain tumors by RNA-Sequencing and ChIP-Sequencing revealed alterations in the molecular landscape associated to H3.3-G34R expression. In particular, the expression of H3.3-G34R modifies the histone marks deposited at the regulatory elements of genes belonging to the JAK/STAT pathway, leading to an increased activation of this pathway. This histone G34R-mediated epigenetic modifications lead to changes in the tumor immune microenvironment of these tumors, towards an immune-permissive phenotype, making these gliomas susceptible to TK/Flt3L immune-stimulatory gene therapy. The application of this therapeutic approach increased median survival of H3.3-G34R tumor bearing animals, while stimulating the development of anti-tumor immune response and immunological memory. Our data suggests that the proposed immune-mediated gene therapy has potential for clinical translation for the treatment of patients harboring H3.3-G34R high grade gliomas.

Developmental exposures to common environmental contaminants, DEHP and lead, alter adult brain and blood hydroxymethylation in mice.

Introduction: The developing epigenome changes rapidly, potentially making it more sensitive to toxicant exposures. DNA modifications, including methylation and hydroxymethylation, are important parts of the epigenome that may be affected by environmental exposures. However, most studies do not differentiate between these two DNA modifications, possibly masking significant effects. Methods: To investigate the relationship between DNA hydroxymethylation and developmental exposure to common contaminants, a collaborative, NIEHS-sponsored consortium, TaRGET II, initiated longitudinal mouse studies of developmental exposure to human-relevant levels of the phthalate plasticizer di(2-ethylhexyl) phthalate (DEHP), and the metal lead (Pb). Exposures to 25 mg DEHP/kg of food (approximately 5 mg DEHP/kg body weight) or 32 ppm Pb-acetate in drinking water were administered to nulliparous adult female mice. Exposure began 2 weeks before breeding and continued throughout pregnancy and lactation, until offspring were 21 days old. At 5 months, perinatally exposed offspring blood and cortex tissue were collected, for a total of 25 male mice and 17 female mice (n = 5-7 per tissue and exposure). DNA was extracted and hydroxymethylation was measured using hydroxymethylated DNA immunoprecipitation sequencing (hMeDIP-seq). Differential peak and pathway analysis was conducted comparing across exposure groups, tissue types, and animal sex, using an FDR cutoff of 0.15. Results: DEHP-exposed females had two genomic regions with lower hydroxymethylation in blood and no differences in cortex hydroxymethylation. For DEHP-exposed males, ten regions in blood (six higher and four lower) and 246 regions (242 higher and four lower) and four pathways in cortex were identified. Pb-exposed females had no statistically significant differences in blood or cortex hydroxymethylation compared to controls. Pb-exposed males, however, had 385 regions (all higher) and six pathways altered in cortex, but no differential hydroxymethylation was identified in blood. Discussion: Overall, perinatal exposure to human-relevant levels of two common toxicants showed differences in adult DNA hydroxymethylation that was specific to sex, exposure type, and tissue, but male cortex was most susceptible to hydroxymethylation differences by exposure. Future assessments should focus on understanding if these findings indicate potential biomarkers of exposure or are related to functional long-term health effects.

Oral Microbiome Community Composition in Head and Neck Squamous Cell Carcinoma.

The impact of the oral microbiome on head and neck cancer pathogenesis and outcomes requires further study. 16s rRNA was isolated and amplified from pre-treatment oral wash samples for 52 cases and 102 controls. The sequences were binned into operational taxonomic units (OTUs) at the genus level. Diversity metrics and significant associations between OTUs and case status were assessed. The samples were binned into community types using Dirichlet multinomial models, and survival outcomes were assessed by community type. Twelve OTUs from the phyla Firmicutes, Proteobacteria, and Acinetobacter were found to differ significantly between the cases and the controls. Beta-diversity was significantly higher between the cases than between the controls (p < 0.01). Two community types were identified based on the predominant sets of OTUs within our study population. The community type with a higher abundance of periodontitis-associated bacteria was more likely to be present in the cases (p < 0.01), in older patients (p < 0.01), and in smokers (p < 0.01). Significant differences between the cases and the controls in community type, beta-diversity, and OTUs indicate that the oral microbiome may play a role in HNSCC.

IndepthPathway: an integrated tool for in-depth pathway enrichment analysis based on single-cell sequencing data.

MOTIVATION: Single-cell sequencing enables exploring the pathways and processes of cells, and cell populations. However, there is a paucity of pathway enrichment methods designed to tolerate the high noise and low gene coverage of this technology. When gene expression data are noisy and signals are sparse, testing pathway enrichment based on the genes expression may not yield statistically significant results, which is particularly problematic when detecting the pathways enriched in less abundant cells that are vulnerable to disturbances. RESULTS: In this project, we developed a Weighted Concept Signature Enrichment Analysis specialized for pathway enrichment analysis from single-cell transcriptomics (scRNA-seq). Weighted Concept Signature Enrichment Analysis took a broader approach for assessing the functional relations of pathway gene sets to differentially expressed genes, and leverage the cumulative signature of molecular concepts characteristic of the highly differentially expressed genes, which we termed as the universal concept signature, to tolerate the high noise and low coverage of this technology. We then incorporated Weighted Concept Signature Enrichment Analysis into an R package called "IndepthPathway" for biologists to broadly leverage this method for pathway analysis based on bulk and single-cell sequencing data. Through simulating technical variability and dropouts in gene expression characteristic of scRNA-seq as well as benchmarking on a real dataset of matched single-cell and bulk RNAseq data, we demonstrate that IndepthPathway presents outstanding stability and depth in pathway enrichment results under stochasticity of the data, thus will substantially improve the scientific rigor of the pathway analysis for single-cell sequencing data. AVAILABILITY AND IMPLEMENTATION: The IndepthPathway R package is available through: https://github.com/wangxlab/IndepthPathway.

Perinatal Lead Exposure Promotes Sex-Specific Epigenetic Programming of Disease-Relevant Pathways in Mouse Heart.

Environmental contaminants such as the metal lead (Pb) are associated with cardiovascular disease, but the underlying molecular mechanisms are poorly understood. In particular, little is known about how exposure to Pb during early development impacts the cardiac epigenome at any point across the life course and potential differences between sexes. In a mouse model of human-relevant perinatal exposures, we utilized RNA-seq and Enhanced Reduced Representation Bisulfite Sequencing (ERRBS) to investigate the effects of Pb exposure during gestation and lactation on gene expression and DNA methylation, respectively, in the hearts of male and female mice at weaning. For ERRBS, we identified differentially methylated CpGs (DMCs) or differentially methylated 1000 bp regions (DMRs) based on a minimum absolute change in methylation of 10% and an FDR < 0.05. For gene expression data, an FDR < 0.05 was considered significant. No individual genes met the FDR cutoff for gene expression; however, we found that Pb exposure leads to significant changes in the expression of gene pathways relevant to cardiovascular development and disease. We further found that Pb promotes sex-specific changes in DNA methylation at hundreds of gene loci (280 DMCs and 99 DMRs in males, 189 DMCs and 121 DMRs in females), and pathway analysis revealed that these CpGs and regions collectively function in embryonic development. In males, differential methylation also occurred at genes related to immune function and metabolism. We then investigated whether genes exhibiting differential methylation at weaning were also differentially methylated in hearts from a cohort of Pb-exposed mice at adulthood. We found that a single gene, Galnt2, showed differential methylation in both sexes and time points. In a human cohort investigating the influence of prenatal Pb exposure on the epigenome, we also observed an inverse association between first trimester Pb concentrations and adolescent blood leukocyte DNA methylation at a locus in GALNT2, suggesting that this gene may represent a biomarker of Pb exposure across species. Together, these data, across two time points in mice and in a human birth cohort study, collectively demonstrate that Pb exposure promotes sex-specific programming of the cardiac epigenome, and provide potential mechanistic insight into how Pb causes cardiovascular disease.

Zinc Finger MYND-Type Containing 8 (ZMYND8) Is Epigenetically Regulated in Mutant Isocitrate Dehydrogenase 1 (IDH1) Glioma to Promote Radioresistance.

PURPOSE: Mutant isocitrate dehydrogenase 1 (mIDH1) alters the epigenetic regulation of chromatin, leading to a hypermethylation phenotype in adult glioma. This work focuses on identifying gene targets epigenetically dysregulated by mIDH1 to confer therapeutic resistance to ionizing radiation (IR). EXPERIMENTAL DESIGN: We evaluated changes in the transcriptome and epigenome in a radioresistant mIDH1 patient-derived glioma cell culture (GCC) following treatment with an mIDH1-specific inhibitor, AGI-5198. We identified Zinc Finger MYND-Type Containing 8 (ZMYND8) as a potential target of mIDH1 reprogramming. We suppressed ZMYND8 expression by shRNA knockdown and genetic knockout (KO) in mIDH1 glioma cells and then assessed cellular viability to IR. We assessed the sensitivity of mIDH1 GCCS to pharmacologic inhibition of ZMYND8-interacting partners: HDAC, BRD4, and PARP. RESULTS: Inhibition of mIDH1 leads to an upregulation of gene networks involved in replication stress. We found that the expression of ZMYND8, a regulator of DNA damage response, was decreased in three patient-derived mIDH1 GCCs after treatment with AGI-5198. Knockdown of ZMYND8 expression sensitized mIDH1 GCCs to radiotherapy marked by decreased cellular viability. Following IR, mIDH1 glioma cells with ZMYND8 KO exhibit significant phosphorylation of ATM and sustained γH2AX activation. ZMYND8 KO mIDH1 GCCs were further responsive to IR when treated with either BRD4 or HDAC inhibitors. PARP inhibition further enhanced the efficacy of radiotherapy in ZMYND8 KO mIDH1 glioma cells. CONCLUSIONS: These findings indicate the impact of ZMYND8 in the maintenance of genomic integrity and repair of IR-induced DNA damage in mIDH1 glioma. See related commentary by Sachdev et al., p. 1648.

Chemoradiation therapy alters the PD-L1 score in locoregional recurrent squamous cell carcinomas of the head and neck.

PD-L1 testing guides therapeutic decision-making for head and neck squamous cell carcinoma (HNSCC). We sought to understand whether chemoradiation therapy (CRT) influences the PD-L1 combined positive score (CPS) and other biomarkers of response to immunotherapy. PD-L1 expression was assessed using immunohistochemistry, and bulk RNA sequencing was performed on 146 HNSCC patients (65 primary sites, 50 paired local recurrences, and 31 paired regional recurrences). PD-L1 was scored using the CPS of ≥1, ≥20, and ≥50. Overall, 98 %, 54 %, and 17 % of HNSCCs had a CPS ≥1, ≥20, and ≥50, respectively. When using a cut-off of ≥1, CRT did not significantly change CPS at the locoregional recurrent site. However, there were significant changes when using CPS ≥20 or ≥50. The CPS changed for 32 % of patients when using a CPS ≥20 (p < 0.001). When using a CPS ≥50, there was a 20-23 % (p = 0.0058-0.00067) discordance rate at the site of locoregional recurrence. Oral cavity cancers had a significantly higher discordant rate than other primary sites for CPS ≥50, 44 % (8/18, p = 0.0058) and 58 % (7/12, p = 0.00067) discordance at the site of local and regional recurrence, respectively. When evaluating the 18 gene IFN-É£ signature predictive of response to anti-PD-1 blockade, there was a statistically significant increase in the IFN-É£ signature in recurrent larynx cancer (p = 0.02). Our study demonstrates that when using a higher cut-off of CPS ≥20 and ≥50, a repeat biopsy may be warranted after CRT for local and regional recurrent HNSCCs.

Comprehensive enhancer-target gene assignments improve gene set level interpretation of genome-wide regulatory data.

BACKGROUND: Revealing the gene targets of distal regulatory elements is challenging yet critical for interpreting regulome data. Experiment-derived enhancer-gene links are restricted to a small set of enhancers and/or cell types, while the accuracy of genome-wide approaches remains elusive due to the lack of a systematic evaluation. We combined multiple spatial and in silico approaches for defining enhancer locations and linking them to their target genes aggregated across >500 cell types, generating 1860 human genome-wide distal enhancer-to-target gene definitions (EnTDefs). To evaluate performance, we used gene set enrichment (GSE) testing on 87 independent ENCODE ChIP-seq datasets of 34 transcription factors (TFs) and assessed concordance of results with known TF Gene Ontology annotations, and other benchmarks. RESULTS: The top ranked 741 (40%) EnTDefs significantly outperform the common, naïve approach of linking distal regions to the nearest genes, and the top 10 EnTDefs perform well when applied to ChIP-seq data of other cell types. The GSE-based ranking of EnTDefs is highly concordant with ranking based on overlap with curated benchmarks of enhancer-gene interactions. Both our top general EnTDef and cell-type-specific EnTDefs significantly outperform seven independent computational and experiment-based enhancer-gene pair datasets. We show that using our top EnTDefs for GSE with either genome-wide DNA methylation or ATAC-seq data is able to better recapitulate the biological processes changed in gene expression data performed in parallel for the same experiment than our lower-ranked EnTDefs. CONCLUSIONS: Our findings illustrate the power of our approach to provide genome-wide interpretation regardless of cell type.

ATRX loss in glioma results in dysregulation of cell-cycle phase transition and ATM inhibitor radio-sensitization.

ATRX, a chromatin remodeler protein, is recurrently mutated in H3F3A-mutant pediatric glioblastoma (GBM) and isocitrate dehydrogenase (IDH)-mutant grade 2/3 adult glioma. Previous work has shown that ATRX-deficient GBM cells show enhanced sensitivity to irradiation, but the etiology remains unclear. We find that ATRX binds the regulatory elements of cell-cycle phase transition genes in GBM cells, and there is a marked reduction in Checkpoint Kinase 1 (CHEK1) expression with ATRX loss, leading to the early release of G2/M entry after irradiation. ATRX-deficient cells exhibit enhanced activation of master cell-cycle regulator ATM with irradiation. Addition of the ATM inhibitor AZD0156 doubles median survival in mice intracranially implanted with ATRX-deficient GBM cells, which is not seen in ATRX-wild-type controls. This study demonstrates that ATRX-deficient high-grade gliomas (HGGs) display Chk1-mediated dysregulation of cell-cycle phase transitions, which opens a window for therapies targeting this phenotype.

Characterization of the immune response in patients with cancer of the oral cavity after neoadjuvant immunotherapy with the IRX-2 regimen.

OBJECTIVE: IRX-2 is a homologous cell-derived multi-cytokine biologic with multifaceted immune modulatory effects that has been shown to induce increased lymphocyte infiltration into primary tumors in oral cavity carcinoma. Our objective was to characterize tumor immune gene expression and epigenomic changes after neoadjuvant IRX-2 immunotherapy in patients with squamous cell carcinoma of the oral cavity. METHODS: A randomized phase II trial was conducted of the IRX regimen 3 weeks prior to surgery for previously untreated patients with Stage II-IV oral cavity carcinoma. The treatment regimen consisted of low dose (300 mg/m2) cyclophosphamide (day 1) followed by 10 days of regional perilymphatic IRX-2 cytokine injections and daily oral indomethacin, zinc and omeprazole (Regimen 1) compared to the identical regimen without the IRX-2 cytokines (Regimen 2). The NanoString immune panel (730 genes) and Infinium MethylationEPIC BeadChip were performed to assess the gene expression and DNA methylation signatures, respectively, in pre- and post-immunotherapy tumor samples. RESULTS: A total of 51 and 79 immune-related genes were found upregulated and downregulated, respectively, in the samples from Regimen 1 patients after treatment, while 51 and 56 were found upregulated and downregulated in the samples for Regimen 2. When comparing the changes between the two regimens, we identified 9 genes significantly different, including DMBT1, a potential tumor suppressor, functioning in tumor invasion of head and neck cancer. The exploration of DNA methylation showed slight overall hypermethylation after treatment in both regimens, especially for Regimen 1 immune responders, and methylation-based cell type deconvolution demonstrated high concordance with tumor infiltrating T lymphocyte cell counts. CONCLUSION: While a consistent patient response after treatment was observed, most changes were similar between regimens, indicating a subtle, targeted, or patient-specific effect of IRX-2 cytokines. Change in DMBT1 expression was a unique finding that will require further study to better understand its significance.

Molecular Tumor Subtypes of HPV-Positive Head and Neck Cancers: Biological Characteristics and Implications for Clinical Outcomes.

Until recently, research on the molecular signatures of Human papillomavirus (HPV)-associated head and neck cancers mainly focused on their differences with respect to HPV-negative head and neck squamous cell carcinomas (HNSCCs). However, given the continuing high incidence level of HPV-related HNSCC, the time is ripe to characterize the heterogeneity that exists within these cancers. Here, we review research thus far on HPV-positive HNSCC molecular subtypes, and their relationship with clinical characteristics and HPV integration into the host genome. Different omics data including host transcriptomics and epigenomics, as well as HPV characteristics, can provide complementary viewpoints. Keratinization, mesenchymal differentiation, immune signatures, stromal cells and oxidoreductive processes all play important roles.

Perinatal DEHP exposure induces sex- and tissue-specific DNA methylation changes in both juvenile and adult mice.

Di(2-ethylhexyl) phthalate (DEHP) is a type of phthalate plasticizer found in a variety of consumer products and poses a public health concern due to its metabolic and endocrine disruption activities. Dysregulation of epigenetic modifications, including DNA methylation, has been shown to be an important mechanism for the pathogenic effects of prenatal exposures, including phthalates. In this study, we used an established mouse model to study the effect of perinatal DEHP exposure on the DNA methylation profile in liver (a primary target tissue of DEHP) and blood (a common surrogate tissue) of both juvenile and adult mice. Despite exposure ceasing at 3 weeks of age (PND21), we identified thousands of sex-specific differential DNA methylation events in 5-month old mice, more than identified at PND21, both in blood and liver. Only a small number of these differentially methylated cytosines (DMCs) overlapped between the time points, or between tissues (i.e. liver and blood), indicating blood may not be an appropriate surrogate tissue to estimate the effects of DEHP exposure on liver DNA methylation. We detected sex-specific DMCs common between 3-week and 5-month samples, pointing to specific DNA methylation alterations that are consistent between weanling and adult mice. In summary, this is the first study to assess the genome-wide DNA methylation profiles in liver and blood at two different aged cohorts in response to perinatal DEHP exposure. Our findings cast light on the implications of using surrogate tissue instead of target tissue in human population-based studies and identify epigenetic biomarkers for DEHP exposure.